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Abstract 

The study concerns a special symbolic calculus of interest for signal analysis. This cal- 
culus associates functions on the time-frequency half-plane / > with linear operators 
defined on the positive- frequency signals. Full attention is given to its construction which 
is entirely based on the study of the affine group in a simple and direct way. The correspon- 
dence rule is detailed and the associated Wigner function is given. Formulas expressing the 
basic operation (star-bracket) of the Lie algebra of symbols, which is isomorphic to that 
of the operators, are obtained. In addition, it is shown that the resulting calculus is co- 
variant under a three-parameter group which contains the affine group as subgroup. This 
observation is the starting point of an investigation leading to a whole class of symbolic 
calculi which can be considered as modifications of the original one. 



1 Introduction 

The notion of time-frequency analysis of signals received a theoretical basis when D.Gabor 
|l| proposed to interpret signals by operations very similar to those of quantum mechanics. 
In this approach, the real signals are replaced by their positive-frequency parts (the "an- 
alytic" signals) and the communication characteristics are introduced as the mean values 
of hermitian operators defined on the corresponding Hilbert space. Actually, the analogy 
with quantum mechanics has rapidly suggested to recognize the time and the frequency 
as two conjugate variables and has led to look for a phase space formulation of the theory. 
Various solutions have then been proposed which were essentially adaptations of solutions 
developed precedently in the quantum mechanical context. 

From a mathematical point of view, the question of the phase space formulation of 
a quantum theory appears as a problem of symbolic calculus since each operator of the 
theory has to be replaced by a function (its symbol) of the phase space variables pj. The 
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construction of such a correspondence rule is determined by the class of operators to sym- 
bolize and by subsidiary constraints relative to the physical interpretation of the formalism. 
In signal theory, this last point concerns simply the formal invariance of the operations in 
changes of reference clocks. The corresponding constraint will be automatically ensured if 
the correspondence rule is covariant with respect to relevant representations of the affine 
group. The object of the following developments is to construct a symbolic calculus on 
the time-frequency half-plane along these lines. 

Section 2 gives a study of the affine group and its subgroups, their representations 
in the Hilbert space TL of signals and their actions in the time-frequency half-plane V. 
Decompositions of TL and V into invariant subspaces are carried out in parallel and shown 
to be related in a natural way. The correspondence rule is derived in Section 3 by exploiting 
the completeness of the invariant decompositions just performed. Some properties of this 
rule are given in Section 4 and the associated Wigner function is determined. Section 
5 is devoted to the derivation of the star product, which allows to symbolize products 
of operators and is the core of the symbolic calculus. In Section 6, the star bracket is 
introduced and shown to reduce to an ordinary Poisson bracket in some special cases. As 
an application of this remark, a Hamiltonian flow is defined in phase space. Finally, it 
is shown in Section 7 that the method has possible extensions when working with larger 
groups containing the affine group. 



2 The affine group in signal theory 

The signals to be considered are real- valued functions of time s(t) which are interpreted 
independently of the choice of time origin and units. This means that the description of 
the signal processing operations must not change if time t is replaced by at + b where b is 
real and a positive. The corresponding change on the signal can be readily written down. 
Notice that there is the freedom of rescaling the signal when the units are changed. The 
transformations on the signal that will be considered are written as: 

s(t) ^a r sia-^t-b)) (2.1) 

where r is a real exponent. 

The set of transformations 
composition law expressed by: 

{a,b)(a',b') = {aa',b + ab') (2.2) 

The group structure is the mathematical expression of the equivalence of the reference 
clocks. To stress the fact that only positive dilations are considered, the parameter u such 
that: 

a = e u (2.3) 

will be used instead of a. 



2.1) forms the affine group, with parameters (a, b) and 
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2.1 Representation in the Hilbert space of positive- frequency signals 

The Fourier transform of the signal s(t) will be introduced by: 

!(/) = / S (t)e~ 2 ^ dt (2.4) 



IR 

Since s(t) is a real- valued function, it is characterized by its positive frequency part only. In 
agreement with the usual practice in signal theory, this property will be used to substitute 
to s(t) the so-called analytic signal whose Fourier transform is given by: 

S(f) = 2Y(f)s(f) (2.5) 

where Y(f) is Heaviside's step function. 

The counterpart of transformation ( |2.lD for the analytic signal S(f) is easily obtained. 
Namely, by action of an element (u, b) of the affine group, signal S{f) goes to a new 
function [U(u, b)S](f) defined by: 

U(u, b)S(f) = e {T+l)u e- 2i7Tfb S(e u f), r € M (2.6) 

The operators ~U(u, b) constitute an irreducible representation of the affine group that 
is unitary for the scalar product defined by 

(S,S')= f S(f)S'*(f)f 2r+1 df (2.7) 

The Hilbert space of functions S(f) on the half-line which are square integrable for 
the measure f 2r+1 df will be denoted by TC (the reference to r being implicit). From a 
physical point of view, it may be necessary to distinguish between signals transforming 
according to ( |2.6| ) with different values of r. However, the corresponding representations 
are easily shown to be unitarily equivalent. In fact, it may be shown that the affine group 
has only two inequivalent unitary representations differing by the sign of the imaginary 
exponential in fl2,6j ) Q. 

The infinitesimal operators of representation (|2.6j) are introduced by the relations: 



=6=0 



= fS(f) (2.9) 
The operators (3 and f thus obtained verify the commutation relation: 

M = ~t (2-10) 

The subgroups of the affine group and their representations will now be studied in detail. 
There are two kinds of subgroups consisting either of pure dilations or of pure translations. 
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The translation subgroup is invariant by conjugation. On the contrary, since a dilation 
can be performed from an arbitrary time origin £, there is a whole family of conjugate 
dilation subgroups G^ which consist of the following set of elements: 

G 5 = {(e",£(l-e"))} (2.11) 
The restriction of representation (^^) to G^ is given by 

U 5 (n)5(/) = e( r+1 )" e- 2i7Tf ^-e u )s{e u f) (2.12) 
and the corresponding generator is 

= G9-e/W) (2.13) 



«=o 



2iir du 

where operators f3 and / are given by fl2.Sj ) and ( |2.9[) respectively. 



The operator of representation (2.12) can be written in terms of generator as: 

U € = e - 2inuf3 t (2.14) 

For a given £;, the eigenfunctions of the self-adjoint operator (3^ are defined by the 
equation 

/3^(/) = /^|(/), /?GB (2.15) 
whose solution is, up to a multiplicative constant: 

= f- 2 ™^- 1 e' 2i7TU (2.16) 

By construction, the functions ipg(f) are also eigenfunctions of the unitary operator ( p. 14 ) 
and verify 

Ut(u)4(f) = e- 2i ^ u 4(f) (2.17) 

For a given £, functions ( 2.16j ) constitute an improper basis of the Hilbert space TL. In 
fact, one verifies that there is orthogonality of the functions corresponding to different (3 

(4(f),4,(f))=6(p-f3') (2.18) 

where the scalar product has the form ( |2.7| ). The completeness of this basis is expressed 
by: 

/ ^(f)4*(f')d(3 = r 2r - 1 6(f-f) (2.19) 

For each £, the set of all the projectors on the functions gives a decomposition of 

the Hilbert space TL. From ( 2.17j ) it results that each of these projectors is invariant by 
the transformations of the subgroup G^. 

More generally, the action of the whole affine group on functions V'iC/) is given by: 

V(u,b)4(f) = e- 2i ^4(f) (2.20) 

with £' = e u £ + b. This shows that the set of projectors associated with subgroup G^ is 
transformed as a whole into the set of projectors associated with subgroup G e u^ + f,. In 
other words, this means that the set of all decompositions of TL by the various subgroups 
Cr£ is invariant by action of the whole group. 
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2.2 Action of the group in the time- frequency half-plane 



The physical phase space for signal analysis is the time-frequency half-plane, T = {(t, /)}, 
t real, / > 0. This space can also be obtained as an orbit of the coadjoint representation 
of the affine group whose action in variables (t, f) is given by: 

(u,b): ( t J)^(e u t + b,e- u f) (2.21) 

where the pair (u, b) characterizes an element of the group. 

This action induces a transformation on functions A(t, f) defined on T. This transfor- 
mation has the form: 

(u, b) : Ait, f) — A(e- U (t - 6), e u f) (2.22) 

According to definition fl2.11| ), the restriction of representation ( p. 21 ) to subgroup G^ 
reads: 

( t>/ )_>( e « t + £(l-e u ), e -'7) (2-23) 
and the corresponding transformation of phase space functions is given by: 

Ait, f) — » A(e~ u (t -£)+£, e u f) (2.24) 

The transformation ( 2.23| ) has the following invariant: 

(t-0/ = A /3elR (2.25) 

where (3 is an arbitrary constant. For each £, this relation defines a family of curves in 
phase space labelled by (3. Actually any point in V belongs to one, and only one, of these 
curves and a partition of phase space has thus been achieved for any value of £. 
Uniform distributions on curves are given by the Dirac distributions: 

4 = £)/ -/3) (2.26) 

These distributions are invariant by action of the subgroup and transform under the 
full group as follows: 

in, b): 5i = 6((t - Of ~ 0) — » <5((e"«(t - b) - $e u f - 0) = 5f +h (2.27) 

This result can be compared with ( |2.20| ). In fact, the full group connects the various phase 
space partitions associated with subgroups G^. 



2.3 Correspondence between invariant structures 

The subgroup G^, considered for any given value of £, has allowed to perform decompo- 
sitions of both 7i and V where the resulting partitions have been labelled by parameters 
denoted respectively by (3 and (3. We will now show that a natural connection does exist 
between these two parameters. To develop this point, it is useful to introduce the notion 
of localized signals in a consistent way. 

A signal St (f) £ T~C is said to be localized at time to provided an affine transformation, 
defined by the element (u, b) of the affine group, sends it into a signal localized at (e u to + b). 
This is expressed by the following equation: 

XJ(u,b)S t0 (f) = e^ u e- 2 ^ b S t0 (e u f) = S e u to+b (f) (2.28) 
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whose solution is, up to a multiplicative constant: 



Sto(f) = r 1 " 1 e-^ ft0 (2.29) 

Notice the factor /~ r_1 which is required by dilation invariance. In this way, the whole 
set of localized states St is stable by action of the affine group. 

In the decomposition of Tt relative to the subgroup Gj, the /3-subspace is generated 
by tpp{f) of the form ( 2. 16| ) . Developing the phase of this function in the vicinity of a 



frequency /o leads to the approximate expression: 

4(f) k e-^^/o- 1 ) f-r-l e -**Mf6+t) (2.30) 



A comparison with ( 2.29[ ) shows that, apart from a constant phase, the approximated 



function can be seen as a state localized at: 

to = P/fo + £ (2.31) 



Thus, in a local study of functions ( 2.16Q , it is possible to associate with any frequency / 



a time t given by t = (3/ f + £. This relation defines precisely an orbit of subgroup in 
phase space V (cf.( 2.25| )). It leads naturally to identify the parameters of the connected 



elements by setting: 

~0 = f3 (2.32) 

This identification, which is clearly invariant under the group action, will play a central 
role in the following construction. 

3 Geometric correspondence between symbols and opera- 
tors 

The aim now is to set up a correspondence between operators on the Hilbert space 7i 
and functions (their symbols) on the time-frequency half-plane V. The approach will be 
essentially based on the parallel decompositions of 7i and Y introduced in the preceding 
section. The derived correspondence will be said to be geometric to emphasize the fact 
that it is essentially founded on a study of the group and of the invariants associated with 
its subgroups. 

3.1 Characterization of operators by their diagonal elements in the var- 
ious £-bases 

Consider an operator A defined by a kernel according to: 

[ASK/0 = / MfiJ2)S(f 2 )f! r+1 df 2 (3.1) 



Its diagonal matrix elements on the G^-invariant basis defined by ( 2.1 6| ) are given by: 
In&P) = W>|,A^) (3.2) 

oo roo 

4*(h)A{hj 2 )4u 2 ){hh) 2r+i dhdh 
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Quantities have several properties that will be detailed below. First, for 

any given value of £ and f3, they are invariant by action of the representation on 
the operator A. Moreover, the operator A can be recovered from the values of I-h(^,P) 
(—00 < £ < 00, —00 < (3 < 00) and a simple formula exists for the trace of operators. To 
be able to discriminate between functions /?) associated with distinct operators, we 
will also use the notation /^(A;£,/3) for (|3,2|). 

G^-invariance o/7^(£,/3) 

The matrix elements of the transformed operator A' = U^" 1 AU^ are given by : 

(V^U^AU^I) = (U^,AU^) 

the last equality being a consequence of property ( |2. 17 ). 
Affine covariance 

The action of the full group on Jft(£, (3) is computed in the same way, using property 
( ^20|) . One obtains: 

I w ( A'; £,/?) = I w ( A; e«£ + &,/?) (3.3) 
where A' = U -1 (tt> d)AU(«, b). 

This relation expresses the diagonal elements of A' in the £-basis in terms of the 
diagonal elements of A in the basis corresponding to the subgroup G^, conjugate of 
G^ by the transformation (u, b) of the affine group. 

Reconstruction of the operator A 

Defining the two-dimensional Fourier transform of P) by 

T n (u, v)= ( e- 2l « + ^l n {^ 0) d£dp (3.4) 



Ju 2 

we can compute it from ([T^) and find 

? \- A ( ue(v/2) uei ~ v/2) \ ( u Y r h 

H[U,V) \2 sinh(«/2)' 2sinh(<;/2) J Usinh(u/2)/ 4sinh 2 (t>/2) [U/V) 

(3-5) 

The inversion of this relation can be done directly by changing the variables. The 
kernel of operator A is thus given by: 

A(/i,/ 2 ) = IhUx ~ /2,M/i// 2 ))(/i/2)- r - 1 |/i - h\ (3-6) 

Trace of a product of two operators 

The scalar product of two operators can be defined as 

Tr(ABt) = / ^A(h,h)B*(h,h){hh? r+l dhdh (3-7) 

where is the hermitian conjugate of operator B. 
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If Ift(A;u,v) and I^(B;ti, v) are the functions corresponding respectively to the 
operators A and B, the computation of the trace of AB^ is performed using (3.6) 
and the result follows: 

Tr(AB t )=/ T n (A;u,v)T n (B*;u,v)\u\ dudv (3.8) 

3.2 Characterization of symbols by their line integrals on ^-orbits 

The description of functions A(t, f) defined on phase space T can be carried out in a way 
which parallels that followed for the operators. The starting point is the association of a 
function Iy(^,/3) with each phase space function A(t,f) by the relation: 

J r (£, 13)= I A(t, f)5((t - Of ~ P) dtdf (3.9) 



For a given £, the function defined by ( |3.9D is interpreted as the integral of A(t, f) on the 
various orbits of the subgroup G^. The study of the properties of function Ir(£,, 0) follows 
the same steps as in Section 3.1. For practical reasons, the notation Ir(A;£,,[3) will also 
be used for expression (|3.9|). 



G?£ -invariance of Ip (£, /3) 

The integral in ( |3.9| ) is clearly invariant under the action of G% on A(t, f) defined 
by (|]D- 

Affine covariance 

When A(t, f) is transformed into A'(t, f) = A(e u t+b, a" 1 /), the function I r (A; f , 0) 
becomes: 

I r (A';tl3) = I r (A;e u t + b,f3) (3.10) 

This relation connects the integrals of the function A' on the G^ orbits to the integrals 
of A on the orbits of the subgroup Gg, conjugate of G^ by the transformation (u, b) 
of the affine group. 

Reconstruction of the function A(t, f) 

Relation ( |3.9| ) defines in fact the Radon transform of Ait, f) with respect to arrays 
of hyperbolas parametrized by £ and (3. To invert this transform, we introduce the 
two-dimensional Fourier transform of Ir(£,/3) in a manner analogous to ( |3.4| ) and 
compute 

I r (u,v)= [ e~ 2i7rut A(t, (u/v)) ^ (3.11) 
Jn \v\ 

The inversion of this formula gives: 

A(t,f) = f [ e^Mfv^M dv (3.12) 



ISR 

Scalar product of two symbols 

Consider two functions A(t,f), B(t,f) and their corresponding functions. The 
use of relation ( |3.12| ) and Parseval's formula yield: 

[ A(t,f)B*{tJ) dtdf = [ T r (A;u,v)T r (l3*;u,v)\u\ dudv (3.13) 
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3.3 The geometric correspondence rule 

So far, operators on the Hilbert space TL and functions on the phase space V have been 
characterized respectively by functions i?i(£, (3) or Ir(£,,(3). Moreover, for each value of £, 
both functions I-^(^,/3) and Ir(£,/3) are invariant by action of the subgroup G%. Thus it 
is possible to define the correspondence between operators and symbols by requiring that 

lH&P)=WS,P) (3-14) 

and to denote the common function by /(£,/?). This relation is stable by the affine group 
as seen from properties (|3.3|) and ( 3.10|) . Actually, the consistency with the group ac- 
tion implies only the proportionality of the two members of (|3.14 ). The interest of the 
strict identification which has been adopted is to ensure that the identity operator will be 
symbolized by the function one. 

According to ( |3.5| ) and fl3.11| ), the Fourier transform of relation (3.14) leads to: 



ft ( Ue(V/2) Uei ' V/2) \ ( — V r H 

[U,V) 1 2 sinh(v/2) ' 2 sinh(t;/2) ) V 2 sinh(i;/2) J 4 sinh 2 (v/2) [U/V) 



e- 2i ^ ut A(t,(u/v)) ^ 
K. \v\ 

The equality of the two right hand sides gives the relation between the kernel A(fi,f 2 ) 
and the symbol A(t,f). This is the correspondence rule that can be rewritten either 
as a formula giving the operator kernel in terms of the symbol (Weyl mapping) or as a 
formula giving the symbol in terms of the kernel (Wigner correspondence). The result is 
the following: 

Result 3.1 The operator A corresponding to the function A(t,f), defined on the time- 
frequency half-plane T, is given by its kernel (defined in fcH. as: 

Conversely, the operator A defined by its kernel A(/i,/2), is represented on phase 
space by a symbol A(t, f) according to: 

A{t,f) -J J^e A y 2sinHv/2 y 2sinHv/2) ) [ 2sinHv/2) ) dv 

In addition, the following unitarity property holds: 

J A(t, f)B*(t, f) dtdf = Tr (ABt) (3.17) 

Among the properties of the correspondence, it can be deduced from ( 3,16| ) that the 
identity operator with kernel: 

/(/i,/ 2 ) = (/i/ 2 )- r - 1/2 5(/i-/ 2 ) (3.18) 
has a symbol equal to one, and that operators such that 

A(f 1 ,f 2 )=A*(f 2 J 1 ) (3.19) 
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are in correspondence with real- valued functions A(t, f). 

At this point, it would be desirable to characterize in a precise way classes of operators 
and functions that are put in correspondence by the above rule. Clearly, because of 
property ( 3,17| ), Hilbert-Schmidt operators can be associated with functions in the space 



L 2 (IR + ; f 2r+1 df). In fact, it is possible to characterize other classes of operators by their 
symbols j§-||, as is done in the usual Weyl case with bounded operators H-[Q]. But this 
is outside the scope of the present paper. 

Result ([O]) concerns only functions defined on the positive half-line and their corre- 
sponding operators. Some side comments concerning its possible extension are given in 
the following remark. 



Remark 3.1 A direct extension of result (|3.l| ) can be obtained if a function A(t, f) defined 
on the whole plane is substituted in ( |3.15| ) and if negative values for f\ and fi are allowed. 
However, the result will be a function A(fi, $2) defined only in the first and third quadrant 
of the (f\, f-z) plane. The regions where f\ and f 2 have different signs are definitely 
forbidden in the above symbolic calculus. 

4 Some special cases 

4.1 Symbol of the projector on a /3-subspace of H 

The /3-subspace associated with subgroup G% is characterized by the function ViCf) i n ~ 



troduced in ( 2.1q ). The projector on this subspace is the linear operator defined by: 

S(f) — (4.1) 



where the scalar product has the form (j2.7| ). The corresponding kernel is given by (cf. 

n|(/i,/ 2 ) = ^(/i)^*(/2) (4.2) 

A direct application of formula ( |3.16| ) shows that the symbols associated with ( [4. 2| ) is the 
distribution S((t — £)/ — (3) with support in T localized on the /3-orbit of subgroup G^. 

4.2 AfRne Wigner function as a particular symbol 

Each signal S(f) defines a real linear functional on the set of hermitian operators by: 

< A >= (AS, S) = {S, AS) (4.3) 



where the scalar product is given by ( |2.7| ). 
This functional can also be written as: 



<A>=Tr(n 5 A) (4.4) 



where the trace operation has been defined by (|3.7|) and where lis is the projector on 
signal S(f) with kernel given by: 

U s (f 1 ,f 2 ) = S(f 1 )S*(f 2 ) (4.5) 
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From the reciprocal correspondence ( 3.15 ), ( 3.16 ) and its properties, it results that each 
real linear functional on the set of hermitian operators can be expressed as a real linear 
functional of the real-valued functions on the time-frequency half-plane. The Wigner 
function for the signal S(f) is defined as the function (or distribution) V(t, f) on the 



time- frequency half-plane which allows to rewrite (4.4) in the form: 



< A > = 



V(t,f)A(t,f) dtdf 



(4.6) 



where A(t, f) is the symbol of operator A. From this definition, it is clear that the form of 
the Wigner function is fundamentally dependent of the symbolic calculus which is adopted. 

The equality of the right-hand sides of (^4) and (4.6) gives a relation of type fl3.17|) 
(unitarity property) provided the Wigner function V(t, f) is identified with the geometric 
symbol of the projector Tls- This makes it possible to obtain the relevant Wigner function 
by applying the correspondence rule (3.16) to the kernel ( f4.5p . The result is the affine 
Wigner function H: 



V{t,f) = f 



2r+2 



2inv 



fig 



JR. 



f ve v/2 

2sinh(u/2) 



S* 



fve' v / 2 \ f v 
2sinh(u/2) J V2sinh(w/2) 



2r+2 



dv (4.7) 



This function can be introduced in various ways. In the above development, it arises as a 
consequence of the general correspondence defined by ( [3,15 ),( |3T6| ). 

Expression Q4.7p can also be obtained directly by writing the unitarity relation ( p. 17 ) 
with one symbol equal to V(t, f) and the other one equal to the localized distribution 
S(t — to) 6(f — fo). This leads to the relation: 



P(t ,/o) = Tr(n s A toJo ) 



(4.8) 



where A-t ,fo * s ^ ne operator with symbol S(t — to) S(f — fo). The kernel of this operator 
is given by ( 3.15|) and reads: 



A(/i,/ 2 ) = (/i/ 2 



,-r-i h - h 



-2i7r(/i-/ 2 )to ^ 



h-h 



(4.9) 



ln/i// 2 VW1//2 

Function V(t, f) can then be obtained by ([4.8|). These results are summarized by: 

Result 4.1 The Wigner function of signal S(f) associated with the geometric correspon- 
dence is given by . It can be obtained either as the symbol of the projector II5 on signal 
S(f), or as the mean value of the operator A t() j whose symbol is 5(t — to)5(f — fo). 



The unitarity property ( |3.17 ) gives directly the formulas: 



P(t, f) dtdf = || s 



(4.10) 



J Vi(tj)v 2 (t,f)dtdf = Tr(n 5l n 52 ) 

= \(Si,S 2 )\ 2 



where V\ and V2 stand for the Wigner functions corresponding to S\ and S 2 respectively. 

Many other properties could be mentioned which are related to the interpretation of 
the Wigner function as a time- frequency representation Q , |l(| . 
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4.3 Exponentials of generators and Weyl's formulation of the correspon- 
dence 

We will now show that the symbol of the operator E UQ1 , = e - 2t7T {uof3+v f ) ^ w j iere p an( j 
/ are the generators defined by (|2.8| ) and (|2.9|) , is given by the function £ UoVo (t,f) = 

e -2in(u tf+v f) _ 

It follows from the developments of Section 2.1 that the operator E Uo „ satisfies the 
following eigenvalue equation: 

e -2iTT(u p+v f)^-v /u ^j,^ _ e -2iTr(3u ^-v /u n ^ (4.11) 

The kernel of E nol , is obtained directly from the spectral decomposition: 

E U0V0 (fuf2)= I e~ 2 ^ ^ 0,U \fi)^' V0/U °{f2) dp (4.12) 

J IR, 

and has the explicit expression: 

iW/1,/2) = e (r+1) "° e W {-^A(vo/u )(e Uo - l)}5(f ie u ° - f 2 ) f^ 1 (4.13) 
Substituting this result in ( 3.16| ) , we can write the symbol of E Uo „ as: 



£u 0V0 (t,f) = I e 2 ^e^ +1 ^exp{-2ivr(/(t;o/ U o)(e^-l)A( W ))} 

x \{-v) e (2r+2)(v/2) 5{e u °X{v) - \{-v)) dv (4.14) 

where, by definition 

\{ v ) = v \ (4.15) 

v 1 2sinh(v/2) v ' 

Since 

5(e u ° X(v) - X(-v)) = -L-d(v + u ) (4.16) 

X(u ) 



the result of integration in ( 4. is just: 

£u v (t,f) = e -2-(«o/*+«o/) (417 ) 

Thus, we have established the correspondence 

Remark that the variable (5 = tf rather than t appears naturally. It will sometimes be 
convenient to use it by performing the change of variables according to: 

A(P,f) = [A(t,f)] t=f3/f (4.19) 



The correspondence (4. IS) allows to write the operator A in terms of its symbol A(t, f) 
by the Weyl formula: 

A = J A(u, v) e -n<«P+»f) dudv (4.20) 
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where A(u, v) is the Fourier transform of A(/3, f) given by: 



A{u,v) = / e 2in ^ +vf) A(PJ) dpdf (4.21) 



The correspondence between functions and operators defined by ( |4.20| ) and ( 4.21 ) has 



a form analogous to Weyl's original rule. In fact, it is what L.Cohen calls Weyls'rule for 
operators (3 and / [Q. It can also be obtained as a by-product of Kirillov's analysis [12] 
as shown by A.Unterberger 0] (see also R.Shenoy and T. Parks K~ 



Remark 4.1 The symbols of the self-adjoint operators (3 and / are immediately seen to 
be (5 = tf and / respectively. The time operator, which is especially important in signal 
analysis, can be defined as the operator with symbol t. A direct computation leads to the 
result: 

t = (l/2)(/- 1 /3 + /3/- 1 ) (4.22) 



or: 



Wr+1/2 +^l (4.23) 



2m \ f dj 

Such an operator is not essentially self-adjoint on the Hilbert space L 2 (1R + ; f 2r+1 df) and 
hence has no conventional spectral decomposition. However, it is known from the defor- 
mation quantization theory on phase space pj that such a difficulty can be easily overcome 
and does not affect the physical interpretation. 



5 Star product on the symbols 

Operators on the Hilbert space Ti have been mapped to functions on the time-frequency 
half-plane T by the rule (|3.16| ). Since the algebra of operators is non-abelian, the product 
of operators is mapped into a composition law of functions that cannot be the ordinary 
product. This situation is familiar from Weyl's calculus where the composition law of 
functions is known as Moyal's product and is invariant by Heisenberg's group. Analogous 
products, invariant by different groups, have been constructed by deformation of the usual 
product and are generally referred to as star products (^-products) [§]. In the present case, 
the composition law of symbols corresponding to the product of operators leads directly 
to a star product that is invariant by the affine group. 

In the following, it will be convenient to use another notation for the operator 
exp{— 2m{uj3 + vf)}. In fact, according to ( [4.1l| ) and ( [4.12| ), the action of that operator 
on a signal S(f) can be written as 

exp{-2ivr(u/3 + vf)}S(f) = e (r+1)n exp{-2ivr/(v/n)(e n - l)}S{fe u ) (5.1) 

Comparing this expression with the definition ( |2.6| ) of representation U, we are able to 



13 



rewrite (5.1) as []: 



exp{-2m(«/3 + vf)} S(f) = U ( u, -(e u - 1) ) S(f) 



(5.2) 



The use of representation XJ(u, b) allows to take advantage of the group law and to write 
the composition of operators as: 



U(u, 6)U(u', b') = U(u + u', b + e u b' 



(5.3) 



It shows directly that the product of two operators of type ( |5.1| ) is an operator of the 
same type. With notation ( |5.2| ), the form (4.20) of the correspondence rule ( |3.15| ) between 
functions and operators can be written as: 



/ e u - V 

A(u, v)XJ[u, v ) dudv 

K 2 \ U 



(5.4) 



where A(u, v) is defined in ( ^2l| ). 

Let A and B be operators given in terms of their symbols A(f5, f) and B{(5, /) by rule 
(|5.4|) ,( ^21 ). Their product is equal to: 



AB 

or, using 



/ e Ul - 1 \ / e U2 — 1 1 
A(ui,vi)B(u2,V2)~U ui,vi U u 2 ,v 2 ) du\dv\du2dv2 



( e Ul - I 

A(ui,Vi)B(u2,V2)V Ui + U 2 ,Vl 

R 4 V Ui 



+ e u 'v 2 



e U2 _ I 
U2 



du\dv\du2dv2 



This product can be rewritten in the form ( |5.4| ) of the correspondence rule if variables 
(n, v) are introduced such that: 



U = Ul + U2, v- 



1 



1 



u 



vi- 



+ e u ' v 2 



u 2 



The result is: 



AB = / C(u,v)\j(u, v ) dudv 

Jm 2 \ u 



with C(n, v) defined by: 



(5.5) 
(5.6) 



C(u, v) 



A(ui,vi)B(u2,V2)S(u - Ui - u 2 ) 



(5.7) 



•/ , e ui -l u e U2 -l\u 1 + u 2 ,, 1! , 
x s [ v - [ vi— — + e ^2— ^— ) eUl+U2 _ 1 ) du 1 dv 1 du 2 dv 2 



1 This result can also be obtained directly by noticing that, in the description of the affine group by 
2x2 matrices, the exponentiation of the Lie algebra element defined by (u, v) is equal to: 



exp 



u 


V 




e" 
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It follows from ( |5.6| ) that C(u, v) is the Fourier transform of the symbol of the operator 
AB. 

To obtain an explicit expression of the star product A(/3, f) * B(f3, /), it suffices now 
to invert the Fourier transform of C(u, v) and to express the result in terms of A((3, /) and 
B(0, /). In a first step, we obtain: 

A(PJ)*B(P,f)= (5.8) 
r -?(ui,«i)B(u2, v 2 )e 2m « Ul+U2 ^ + ( Vi+v ^Ve- 2 ™f h( - Ul ' Vl > U2 ' v ^ du 1 dv 1 du 2 dv 2 
where the function h is defined by: 

h{u 1 ,v 1 ,U2,v 2 ) = {v 1 +v 2 ) - — — [v l + e l v 2 (5.9) 

e ui+u 2 — I \ Ul U 2 J 

or, after some manipulations: 

%i, v u u2, v 2 ) = (v lU2 - mm) MlM2 ( e u 1+ « 2 _ i) ( 5 - 10 ) 

The exponential of function h(u\, v\, u 2 , v 2 ) in ( |5.8| ) is what makes the difference with a 
classical product. It can be replaced by its series expansion so that the integrals in (5.8) 
can be performed. In the process, the following replacements take place: 

V; ► df., Ui — > — !— dg., i = 1,2 (5.11) 

The development of the star product can thus be written as: 

A(P, f) * B(f3, f) = A((3, f)B([3, f) (5.12) 

+ Y.-Aj-( d hd^-d h d B2 )T(d Bl ,d B2 )) A((3x,fi)B((3 2 ,f 2 ) | /l=/2=/ 
B=1 n! 7 /9i=A=/3 

where the contribution of operator T will be discussed below. 

If ( 5.12j ) is rewritten in terms of t rather than f3, using the definition: 

p = tf, d 3 = (i/f)d t (5.13) 

the result is reminiscent of the Weyl product. Actually, the difference lies in the presence 
of the operator T(d Bl , dp 2 ) which comes from the development in powers of u\ and u 2 of: 

T( \ n ie ^(e^-l)-n 2 (e^-l) 

T( Ul ,u 2 ) = 2 5.14 

uiu 2 (e Ul+U2 — 1) 

Operator T introduces derivatives of all orders with respect to fi\ and (3 2 . At the lowest 
orders, one has 

T(d /3l ,d 32 )= (5.15) 
1 + J^( d fo - 5 /3i) + I2(27r) 2 ^ lgfe + 720(27r)3 ^ft ~ d ^ d k + d h + 59 Px d fo) + ' ' ' 
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As in Weyl's case, the first term in ( |5.12| ) is the usual product and the following terms 
introduce derivatives of the symbols. It can be observed that the star product of two real- 
valued functions is not real. This is a direct manifestation of the fact that the product of 
hermitian operators is not hermitian. 



Remark 5.1 From Remark p.l| , it follows that in a product of operators, the positive and 
negative frequency parts are multiplied separately. Hence, the same property will hold for 
the star product if its expression is extended to symbols defined on the whole plane. 

The star product simplifies when one of the symbols depends only on one variable. 
This will be illustrated by computing its expression when A(f3, f) is equal either to (3 or 
to a function of / alone and when B(f3, /) is the generic distribution: 

A/WcXA /) = " PoW - /o) (5-16) 



with Fourier transform ( 4.21 ) given by: 

A Mo (u,v) = e-^Mo+^/o) (5 . 17) 

(i) Case where A(/3,f) = (3 

The Fourier transform of (3 is : 

A(u,v) = — —S'Cu)SCv) (5.18) 
2m 



Computation of the star product from ( |5.8| ) and ( 5.10 ) gives: 

/3*<5(0-A)W-/o)= (5-19) 
~m J e- 2 ^M/3o-/?)+* 2 /o] exp{2m[« 1( 5 + fv 2 ^^l]} d Ul du 2 dv 2 

where \{u), defined by ( [4.15 ), is such that: 

\{-u) = e - u X(u) (5.20) 
Performing the ui-integration, we can write d5lSj ) as: 
P*6(J3- (3 )5(f - f ) = [35(f3 - (3 )5(f - /„) (5.21) 



+ / / e ^[MP-M+v 2 U-h)\ fA(in(A(« 2 )))) V2 du 2 dv 2 
Jn 2 \du 2 y "') 



(ii) Case where A(f3,f) depends only on f. 

The Fourier transform of A((3,f) = g(f), with g(f) an arbitrary function, is equal 
to: 

A(u,v) = 6(u) [ g(f)e- 2 ^df (5.22) 



Substitution of the Fourier transforms ( |5.17] ) and ( |5.22j ) in the expression of the star 
product (|5^) -(l5TC|) gives: 

g(f) * S{/3 - o )S(f - / ) = I g(f\(-u))e 2i ^-^ + ^-^ dudv (5.23) 
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6 Star bracket of symbols and Hamiltonian flow 



6.1 Definition of the star bracket 

A Lie bracket is now denned on operators by the operation: 

- 2vri[A, B] = -2vri(AB - BA) (6.24) 

This bracket defines a composition law on the subspace of hermitian operators. When 
transformed by the correspondence rule, using the expression for the star product, this 
bracket gives rise to an operation on real-valued functions that will be called star bracket 
and denoted by { , }*. Thus the star bracket of symbols A(f3,f) and B(f3,f) is defined 
by: 

{A, BUP, f) = -2m(A(P, f) * B(J3, f) - B((3, f) * A((3, /)) (6.25) 
According to Q5.12]) , the bracket can be written as; 

{A(P,f),B(PJ)h = [f{d h d P2 - d f2 d Pl )(l + T^dfrdfy + ■■■) (6.26) 

xA(p 1 ,fi)B(J3 2 j2)\ h =f a = / 
The first term in / is a Poisson bracket given by: 

{A(f3,f),B(P,f)}p-f[—Qp gj g~ f Qp-) (6-27) 

In fact, if the variable (3 is replaced by its definition ( ^.13 ), the expression Q6.27 ) 
becomes identical with the usual Poisson bracket on T. This is the same bracket as 
obtained by Kirillov's theory Jl2| when considering phase space as an orbit of the coadjoint 
representation. 

The form ( 6.26| ) of the star bracket is reminiscent of that obtained with the Weyl 
calculus. However, due to the fact that operator T{dp 1 ^dp 2 ) has no symmetry, derivatives 
of arbitrary order will in general appear in the expression of the bracket. 

Some symbols are such that their star bracket with any other function has the form 
of a Poisson bracket. In deformation quantization, such symbols are called preferred 
observables 0], [15]. A typical example is the harmonic oscillator Hamiltonian on 1R 2 
with Weyl quantization. 

In the present case, the preferred observable are represented by the symbols (3, f and 
In/. Explicitly, their star brackets with any symbol X{j3, f) are given by: 

{A #(/?,/)}* = f^-MPJ) (6.28) 
of 

{/,*(£,/)}* = -fJLxtfJ) (6.29) 

{In/, #(0,/)}* = ~§g X ^f) ( 6 - 3 °) 
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The derivation of these formulas can be based on expressions flpip , O) and on their 
counterpart with reverse order in the star products. Final results follow from the relations: 



A(w) - X(-u) = u, X(-u) = e~ u \{u) (6.31) 
which are verified by the function ( [4.15| ) 

6.2 Extended covariance and Hamiltonian flows 

Consider a one-parameter group of transformations on signals S(f) defined by an operator 
O such that: 

S(/; a ) = e- 2 ™°S(/;0) (6.32) 
with S(f;0) = So(f) given. Signal S(f;a) verifies the equation: 

dS ^' 7 = -2ma O S(f; a) (6.33) 
da 

and the projector II on S(f; a) verifies: 

dll 

— = -2ivr[0,n] (6.34) 
oa 

Now, going to phase space, we know that the geometric symbol of II is the affine Wigner 
function (|4.7| ) for signal S(f;a), and that the symbol of the bracket — 2iir[ , ] is the star 



bracket { , Thus equation ( |6.34 ) becomes the following equation in T: 



dV{t d f a ,a) = iO(t, f) , V(t, /; «)K (6-35) 

where 0(t, f) is the symbol of operator O. 

Consider now the case where the symbol 0(t, /) has the particular form: 

H (t,f) = /jtf + vf + alnf (6.36) 



where /i, v and a are real constants. From the values of the star brackets ( 6.28[ ), ( |6.29| ) 



and ( |6.3Q[ ) , it follows that the bracket in ( |6.35| ) reduces to Poisson's bracket. As a result, 
equation ( 6.35| ) takes the form of a Liouville equation with Hamiltonian equal to On(t, /)'■ 

(rPiLJ:a) {^f + uf + alnf,V(t,f;a)} P (6.37) 



da 



(if [Jit + v+ (cr/f))- 



df ^ y IJ " dt 

This means that the a-evolution of the Wigner function is identical to that of an incom- 



pressible fluid in the time- frequency half-plane. Integration of equation ( 6.37 ) with the 
initial condition 

V(tJ;0) = V (t,f) (6.38) 

gives: 

V(t, /; a) = V ((-vfi- 1 + (t + vyT x ) e~" a - af' 1 a) , / e» a ) (6.39) 
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The evolution of the Wigner function in terms of parameter a has a counterpart in Hilbert 
space that can be obtained using the symbolic calculus. The operator Oh for this evolution 
is found directly from (4.20)-( 4.21[ ) and has the expression: 



O h = ///3 + i/f+o-lnf 



(6.40) 



where the operators (3 and f have been introduced in (|2.8| ) and ( |2.9[) . The expression 
( 6.32[ ) written for the operator Oh then gives the a-evolution in the Hilbert space which 
corresponds to the phase space evolution ( |6.39| ). In fact, this leads to a special form of a 
general transformation which can be written as: 



(u,b,c):S(f)^S'(f) 



e (r+l)u e -2iirbf 



f-^ c S(e u f), u,b,c£M 



(6.41) 



This transformation is a projective unitary representation of a three-parameter group Go 
which has been studied in [16]. It has been proved that transformations ( 6.41[ ) correspond 
to symplectic transformations in the phase space. This extended covariance of the calculus 
is at the basis of a generalization of the construction process that will now be considered. 



7 Other symbolizations with geometrical features. 

The method presented above to construct a symbolic calculus on the affine group can be 
embedded in a more general approach based on the systematic consideration of groups 
containing the affine group and having the time-frequency half-plane as coadjoint orbit. 
Such groups belong to the class of simply connected Lie groups having a coadjoint orbit of 
dimension smaller or equal to two [18|. In this class, the groups of interest are »SL(2,IR) 



and the three-parameter solvable groups with Lie algebras defined by generators X a , X b 
and X c satisfying |L7| : 

[x a ,x b ] = x b , [x a ,x c } = kx c , ken (7.42) 

[X a , X b ] = X b , [X a , X c ] = X C + X b (7.43) 

Examples of symbolic calculi invariant by SX(2,]R) have been thoroughly studied in [p^ ]- 
[20]. Here, we will treat the cases of groups Gk with Lie algebras defined by relations 



( 7.42 ) and ( |7.43|) . They provide a natural generalization of the group Go considered so 



far in this paper and will be associated with new forms of symbolic calculus. All groups 



Gk have been introduced in [16] as a tool to find explicit forms of time- frequency pseudo- 
distributions. We refer to that paper for the details concerning the precise construction 
of their unitary representations (projective or not) according to the values of k. Here we 
need only recall the results corresponding to values of k different from and 1. 
The groups Gk, k ^ 0, 1 consist of elements (u,b,c) with the composition law: 

gcf = (u + u',b + e u b f ,c + e ku c') (7.44) 



Their unitary representations in the Hilbert space TL defined in Section 2.1 are given by: 
U k {u, b, c)S(f) = e( r+1 )" e - 2i ^ b}+c ^ S{e u f) (7.45) 



The notation G\ will designate the group with Lie algebra defined in ([7.4 
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The action of the groups in the time-frequency half-plane T is: 

(t, /) — ► (e u t + b + kce- {k - 1)u Z*" 1 , e" u /) (7.46) 

It is clear from ( [7.45 ) and ( |7.46j ) that the cases k = and A; = 1 are degenerate and 
deserve a special treatment (see 

Following the same pattern as in section 2, we determine the subgroups of that 
are conjugate to the dilation subgroup. They are labelled by two real parameters £, t? and 
consist of: 

G^ = {(e u ,e(l-e u ),7 ? (l-e fcM )} (7.47) 
The restriction of representations (|7.45| ) to subgroups G^ v is: 

XJtoS(J) = e Kl)« e -2»(e(i-e»)^(i- e k »)/*) S(/e «) (7 . 48) 
The eigenfunctions of are: 

^(f) = /-2for/9-r-l (7.49) 

They transform under the full group Gk as: 

U(«, 6, e)</f (/) = e- 2 ^" ^ e " +Mefe " +c (/) (7.50) 

In phase space T, the curves invariant by action ( 7.46| ) restricted to G^ v are given by: 

(t-0f-kr]f k = P (7.51) 

A local study of ipp(f) analogous to that performed in Section 2^ leads to the identifica- 
tion: 

= /3 (7.52) 

As in Section 2, we have achieved a correspondence between invariant structures in 7i 
and T. But we now have one extra parameter labelling the vectors ippif) and the curves 
(7.51). To proceed with the construction, we must fix one of these parameters in such a 
way that the affine covariance of the procedure is preserved. Because of the transformation 
law (fr.50D , the only possible choices are either 77 = or f3 = 0o fixed. The case rj = 
yields the calculus just developed in the preceding sections. The only possibility to obtain 
a new calculus is to fix the value of (5 and to let 77 free. 

The same geometric construction as in Section 3 will be carried out. Consider an 
operator A on Ji and form its diagonal matrix elements on V , ^|J(/) : 

Consider a function A(t,f) and form its integral with respect to curves ( [7.51|) : 

I?°(Z,v)=J r A(tJ)5((t-Z)f-krif k -P )dtdf (7.54) 

The study of the properties of functions 1^ (£, 77) and ip° (£, 77) follow the same steps as in 
Section 3. 
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G^-invariance of l£,° (^,tj) and ip (£, rj). 

It can be readily verified, using ( |747D , (|OC| ) and the definitions of 77) and 



Gfc- covariance 

If A' = U^T 1 (n, 6, c)AU/j(m, 6, c) is the operator transformed from A by the group 
representation ( |7.45| ), it follows from ( 7.50|) that: 



J# (A'; £, V) = 1% (A; e u i + b, e ku r, + c) (7.55) 

On the other hand, if A'(t, f) = A(e u t + b + kce~( k ~^ u f k ~ x , a -1 /), a direct com- 
putation leads to: 

1^ (A'; e, r?) = 4° {A; e u £ + b, e ku r? + c) (7.56) 

Possibility of reconstruction of the operator A and of the function A(t, /). 

In spite of the fact that the functions ya Q (f) do not form a basis, it can be shown 
that the operator A is completely characterized by its diagonal elements (|7.53|). 



The reconstruction of the function A(t, f) could be performed by inverting the Radon 
transform ( [7.54] ). 

These properties allow to base the correspondence rule on the identification: 

/# (£, r,) =Ip(Z, V ) (7.57) 

which can be seen as an extension of fl3.14|) . The resulting correspondence rule has the 
form: 

A(t,f)= (7.58) 

f 2r+2 f e 2inup e 2wr(t/-/J )(A fc («)-A fc (-t,)) A(f\ k (u), f X k (-u)) (X k ( U ) X k (- U )) r+1 du 
J IR 

where the functions X k (u) are defined by 

X t (u) = e^(k S ^l) mk " ) (7.59) 



sinh(/cu/2) , 

Conversely the kernel of the operator A can be expressed in terms of the symbol by: 

A(h,f 2 ) = j 'e-^Mjl/i-A-r 1 / 14 ^-^)) (7.60) 

/ f \ -2i7r/3 

x £ (fif2)- r r k \ft 1 -fi' 1 \Atj)dtdf 



J2. 

The operator A corresponding to this kernel can be written in Weyl's form as: 



A= f A(k,p ;u,v)e- 2i ^ uf3+vf ^ dudv (7.61) 
Jtr 2 
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where 



A(k, (3q; u, v) 



,(k-l)u 



(1 - k)u 



fc+l e 2i7r/3o(u-A fc («)+A fc (-u)) 



(7.62) 



e 2iirtf(\ k (u)-X k (-u)) e 2i7Tf(v/u)X k (-«) (e* 



~ l) A{tJ) fdtdf 



It may be noticed that the limit for k = of function Xk(u) is identical to X(u) defined in 
( 4.15| ). In that case, the relation ( [7.62 ) becomes a Fourier transform with respect to /? = tf 
and / and the rule ( 4,20 )-( f4,2l|) is recovered. In the general case, the Fourier transform 
( 4.21D is replaced by the fc-dependent transform ( [7.621) for which Parseval's formula does 
not hold. As a consequence, the unitarity property ( |3.17| ) cannot be obtained. 

However, it is possible to recover an expression for the scalar product of operators in 
terms of symbols. This is accomplished by associating two symbols with any operator A: 
the symbol A(t, /) defined by ( 7.58] ) and a dual symbol Ad{t, /) defined by: 



MtJ) 



(7.63) 



2r+2 



IB 



JKwfh e 2m(t/-/3 )(A fe («)-A fe (- U )) A(f X k (u) , f X k (-u)) fl d (u) du 



where 



d 



du 



to(u) = (X k (u)X k (- U )) r+1 
It can be observed that it is only in the case k = that: 

^(Afc(n) - X k (-u)) = 1 

leading to 



(Xk(u) - X k (-u)) 



(7.64) 

(7.65) 
(7.66) 

(7.67) 



where A(t, f) is the symbol of A and Bd(t, f) is the dual symbol of B. 

These different correspondence rules may be of interest for the study of special opera- 
tors. The case k = —1 where 



A(tJ)=A d (tJ) 
For an arbitrary value of k, the unitarity formula is replaced by: 

Tr(ABt)= / A(t,f)B* d (t,f)dtdf 



A-i(u) 



,u/2 



(7.68) 



stands out. The symbolic calculus corresponding to k = — 1 and /3q = has been intro- 
duced directly by A.Unterberger who applied it in a mathematical context. 

A Wigner function V^°{t,f) can be associated with each fc-calculus. Because of the 
form of relation ( 7.67 ), the Wigner function of signal S(f) is defined as the dual symbol 
of the projector on the signal. We thus associate a time-frequency distribution with each 
value of (k,(3o). When (3q = 0, these distributions coincide with the generalized passive 
distributions introduced in |f[q 1. The case k = 1, which has not been treated explicitly 
here, can be deduced by continuity. 
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8 Conclusion 



A correspondence rule between functions on the time-frequency half-plane and operators 
on the Hilbert space of positive-frequency signals has been obtained. The rule, which is 
said geometric, is entirely based on the study of the affine group and more specially on 
the representations of its subgroups in the two domains. It is obtained, relatively to each 
subgroup, by identification of the decomposition in invariant subspaces with the corre- 
sponding tomographic decomposition in the time-frequency half-plane. The construction 
is very natural and the result can be considered as an affine version of the Weyl rule 
introduced in quantum mechanics. 

In the correspondence, a Lie algebra of hermitian operators is transformed into a Lie 
algebra of symbols, also called star algebra. The basic operations of this algebra are the 
star product and the star bracket of symbols. Formulas for these operations have been 
given. A special attention has been paid to the case where the star bracket reduces to a 
Poisson bracket. 

The obtained symbolic calculus has been integrated into a family of calculi covariant 
by three-parameter groups Gk, k £ 1R, containing the affine group. In this family, the 
geometric rule corresponds to the special value k = and is the only one to ensure the 
unitarity property. 
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